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The ratchet effect is demonstrated theoretically for the simple model of a vortex in a thin superconducting film 
interacting with a periodic array of magnetic dipoles placed in the vicinity of the film . The pinning potential 
for the vortex is calculated in the London limit and found to break spatial inversion symmetry and to depend 
on the orientation of the magnetic dipole moments. The motion of the vortex at zero temperature driven by a 
force oscillating periodically in time is investigated numerically. Drift vortex motion consisting of displacements 
by a translation vector of the dipole array during each period of oscillation is obtained and studied in detail. 
The direction of drift differs in general from that of the driving force, except if the driving force oscillates in 
a direction of high symmetry of the dipole array. The vortex drift velocity depends on the orientation of the 
magnetic moments, and can be tuned by rotating the dipoles. It is pointed out that if the magnetic moments 
are free to rotate, the ratchet effect can be produced and tuned by a magnetic field applied parallel to the film 
surfaces. 



1. Introduction 

Possible applications of superconductivity to 
electronic devices based on the control of vor- 
tex motion has received a great deal of attention 
lately. One line of work uses vortex pinning by 
periodic potentials lacking spatial inversion sym- 
metry to produce drift vortex motion in a prefered 
direction when driven by a oscillating force with 
zero average. This is the so called ratchet effect 
P . Applications to removal of trapped flux from 
superconductors, and to voltage rectification have 
been proposed 2 3 4 5 6 7. Experimental realiza- 
tions of the ratchet effect for vortices have been 
reported by some groups ' 8l9ll0j . The work so 
far has concentrated on vortex pinning potentials 
resulting from periodic modulations of the film 
thickness or from periodic arrays of anti-dots or 
blind holes with asymmetric shape. The objective 
of this paper is to show that the ratchet effect can 
be obtained for vortices in thin superconducting 
films pinned by periodic arrays of magnetic dots 
placed in the vicinity of the film. Here a simple 
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model is studied in detail. One vortex in a thin 
superconducting film interacting with a periodic 
array of point magnetic dipoles, placed on top 
of the film in the London limit. The ratchet ef- 
fect results from the interaction between the vor- 
tex and the magnetic dipole array if the magnetic 
moments are parallel to the film surfaces, and all 
moments point in the same direction. As shown 
here, the interaction between the vortex and the 
dipole array breaks inversion symmetry, and de- 
pends on the orientation of the magnetic moment 
with respect to the lattice. The latter allows the 
ratchet effect to be tuned by rotating the mag- 
netic moments. If the dipole array is made of 
freely rotating magnetic moments, the ratchet ef- 
fect can be produced and tuned by a magnetic 
field applied parallel to the film surfaces. This 
field orients all magnetic moments in the same 
direction, and does not influence the vortices be- 
cause the film is thin. This may be of practi- 
cal interest because it is possible to fabricate ar- 
rays with freely rotating magnetic moments, as 
demonstrated recently by Cowburn, et. al. 11 . 
These authors reported on the magnetic proper- 
ties of arrays of nanomagnets made of Supermal- 
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loy, each nanomagnet being a thin circular disk of 
radius i?, and found that for R ^ 50 — lOOnm the 
magnetic state of each nanomagnet is a single do- 
main one, with the magnetization parallel to the 
disk plane, and that the magnetization can be re- 
oriented by small applied fields. They concluded 
that each nanomagnet acts like a giant magnetic 
moment free to rotate. 

The calculations carried out in this paper start 
from the exact interaction potential in the Lon- 
don limit between a vortex in a thin supercon- 
ducting film and the dipole array. The interac- 
tion potential depends on the orientation of the 
magnetic moments with respect to the dipole lat- 
tice and, in general, lacks spatial inversion sym- 
metry. The equation of motion for the vortex 
driven by a force periodic in time is then solved 
numerically at zero temperature, both for sinus- 
wave and square- wave time dependencies. It is 
found that drift vortex motion takes place in such 
a way that during each period of oscillation the 
vortex displacement is equal to a translation vec- 
tor of the dipole array. The vortex displacement 
is not, in general, in the direction of the driving 
force, and depends in a complicated way on the 
orientation of the magnetic moments and of the 
driving force with respect to the dipole lattice. 
The ratchet effects is, thus, two-dimensional in 
general. One-dimensional ratchet effects are ob- 
tained if the driving force oscillates in a direction 
of high symmetry of the dipole array lattice. In 
this case the vortex motion is found to be simi- 
lar to that reported in the literature for particles 
interacting with one-dimensional ratchet poten- 
tials Q]. The vortex drifts in the direction of the 
driving force oscillation. Its displacement dur- 
ing each period of oscillation is a integer multiple 
of the dipole lattice period in this direction, and 
depends on the orientation of the magnetic mo- 
ments with respect to the direction of the driving 
force. The paper argues that the results obtained 
for this simple model are applicable to thin low- 
Tc superconducting films, and can be extended 
to vortex lattices commensurate with the dipole 
array. 

The paper is organized as follows. In Sec. [2l 
the interaction between the vortex and the dipole 
array is obtained. The motion of the vortex inter- 
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Figure 1. a) Schematic side view of the supercon- 
ducting film with one vortex at r and the dipole 
array on top. b) Top view of the dipole array, 
and definitions of angles a and d 



acting with the dipole array is examined in SecO 
where the main results of the paper are reported. 
Finally, Sec. ^ discusses the limits of validity of 
the model and its possible extensions, and states 
the conclusions of the paper. 

2. Vortex- dipole array interaction 

The superconductor film is assumed to be pla- 
nar, with surfaces parallel to each other and to 
the X — y plane, isotropic, characterized by the 
penetration depth A, and of thickness d <C A. A 
plane dipole array is located above the film at a 
height 2^0 > 0. The dipole array is characterized 
by a rectangular unit cell, with sides ax and a^y, 
and one dipole per cell. All cells have the same 
magnetic moment, m, parallel to the film sur- 
faces, and oriented at an angle a with respect to 
the 3y-axis, as shown in Fig.^ In the case of freely 
rotating magnetic moments, a magnetic filed H 
must be applied parallel to the film surfaces, at 
an angle a with the x-axis, in order to orient all 
magnetic moments in the same direction. The 
magnitude of the field must be large enough to 
avoid reorientation of the dipoles by the fields cre- 
ated by the vortex, by the other dipoles, and by 
the screening and transport currents fiowing in 
the film The field H does not infiuence the 
vortex, because the film is thin <C A). 
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The London limit interactions between vortices 
in superconducting films and magnetic dipoles lo- 
cated outside the film have been discussed in de- 
tail in Ref. 12 , based on the exact solutions of 
London equations obtained in Ref. 13 . The in- 
teractions result from the action of the screening 
current induced by the inhomogeneous field gen- 
erated by the dipoles on the vortices, and can be 
expressed as the magnetostatic interaction of the 
dipoles with the magnetic field generated by vor- 
tices. According to these results, the energy of 
interaction for a single vortex, with unity vortic- 
ity, located at position r, and the dipole array, 
can be written as 

U,-da{r) = - ^ m • bl(R, - r) (1) 

where N is the number of dipoles in the array, Hj 
their positions and is the component parallel 
to the film surfaces of the field generated by the 
vortex. The vortex field is given in the thin film 
limit and for r <C A = 2X^/d by ^ 

(1)0 d r zq 
Mr) = -^^(l-^^^). (2) 

For r > A , b^(r) ~ r/r^. For a periodic dipole 
array, the most important contribution to Uy-da 
comes from dipoles close to the vortex. The in- 
teraction between the vortex and a dipole located 
at distances larger than A from it behaves as 
(R^- - r)/ I R^- - r p, thus falling off as l/i?| for 
Rj ^ r. However, the singular ^/R^ dependence 
is canceled out in the sum over dipoles present in 
Eq. Q by symmetry. The leading contribution 
falls only as 1/R^. Here the contribution from far 
away dipoles is cutoff exponentially, as discussed 
shortly. 

The vortex pinning potential can be written as 
the sum of contributions from each dipole, that is 

Uy-da{r) = - ^i^j-^). (3) 

where U{r) is the pinning potential for a single 
dipole located at the origin, and given by 




Figure 2. Vortex pinning potential (in units of 
Cod = (0o/47rA)^(i) for a single dipole with m par- 
allel to the X-axis, x and y in units of ^. Param- 
eters: m — 0.25^0^0, d = zo = 2(_^ X = 10^. 



where the exponential factor is the cutoff for the 
contributions from far away dipoles. According 
to Eq. (@|), Z//(r) is anti- symmetric with respect to 
spatial inversion and has a minimum (maximum) 
located at r • m = — (+)1.3zo. For m oriented 
parallel to the x-axis, U{r) has the spatial depen- 
dence shown in Fig. [21 

For a dipole array Uv-da{'^) is the periodic rep- 
etition of Z//, centered at each dipole. The spatial 
dependence of Uy-da{'^) for a typical rectangular 
dipole array is shown in Fig. [2l for some values of 
a. 

The vortex-dipole array interaction Uy-da{'^) 
lacks spatial inversion symmetry and leads to the 
ratchet effect when the vortex is driven by an os- 
cillating force, as discussed next. 

3. Ratchet effects 

The motion of a vortex in the potential 
Uy-da{'^)^ driven by a force FL{t) is governed, at 
zero temperature, by 

dr 

rj-=FUt)-VUy-da{r) , (5) 

where rj is the vortex friction coefficient. The 
driving force, FL(t), results from a transport cur- 
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Figure 3. Vortex pinning potentials (in units of 
Cod) for dipole arrays with magnetic moments 
oriented at angle a with the x-axis. x and y 
in units of ^ =vortex core radius. Parameters: 
m = 0.25^0^0, d = zo = 2^, X = 10^. 



rent density, J{t) applied to the film, that is 
Fl(0 = (0o<^/c)J(t) X z. The dipole array lat- 
tice is assumed rectangular with ax > ay ^ zq. 
In this case Uy-da{'^) is like that shown in Fig. 
01 It is convenient to use the following natural 
units for physical quantities. Length: ^ = vor- 
tex core radius. Force: eo = (^o/47rA)^. Time: 
r = r](_/eo. Velocity: cq/t]. Current density: 
Jd = c0o/(12>/37r^A^^), where Jd is the depair- 
ing current. When Eq. Q is written in terms of 
these natural units it depends only on the scaled 
parameters m/(/)o^o, d/zo^ A/zq, J/Jd, and on 
(Rj — r)/zo. The values of J are, of course, lim- 



ited to J < Jd, but in the results reported next 
regions where Jc > Jd are discussed for the sake 
of completeness. 

As will be seen shortly, the dc critical current 
to depin a vortex from a minimum of Uy-dai'^) 
plays an important role in the ratchet effect. To 
obtain it, Eq. is solved for a time-independent 
transport current. The critical current depends, 
in general, on the orientations of m and of the 
driving force with respect to the dipole lattice. 
In the case of interest here, ax, ay ^ zq, it is 
found that the critical current is essentially iden- 
tical to that for a vortex pinned by a single dipole, 
obtained in Ref. 14 , and depends only on the 
angle between the dipole and the driving force, 
/3. The critical current Jc{/3) is shown in Fig. ^ 
The maximum Jc occurs when m is parallel to the 
direction of the driving force, = 0, and the min- 
imum for P = 180^. There is a tenfold difference 
between the maximum and minimum values, and 
a smooth decrease in Jc with increasing /3. The 
maximum and minimum Jc can be estimated ana- 
lytically from the single dipole pinning potential, 
Eq. The result is Jc/Jd - 4m/0o^o for (3 = 0, 
and Jc/Jd - OAm/(t)oZo for p = 180^. 
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Figure 4. Critical current, Jc, for a vortex pinned 
by the dipole array vs. angle between the mag- 
netic moments and the driving force, P (inset). 
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For an oscillating transport current with zero 
average, Eq. Q is solved numerically and the vor- 
tex velocity averaged in time in the steady-state 
regime, denoted by V, is calculated. Drift vortex 
motion occurs if V 7^ 0. The transport current 
is assumed to oscillate in a direction defined by 
the unit vector n, that makes an angle 6 with the 
?/-axis, and to be periodic in time with period P 
and zero average. Writing J(t) = JTf{t)h^ Jt is 
the amplitude, and fit) is chosen as either as si- 
nus wave, f{t) = sin(27rt/P), or as square wave, 
f{t) = 1 for kP <t < {k^ 1/2)^, and f{t) = -1 
for {k - 1/2)P <t < kP, with k = 0, ±1, ±2, .... 
The driving force oscillates in the direction per- 
pendicular to n, which makes an angle with the 
X-axis (Fig.^. It is found that for arbitrary the 
ratchet effect is two-dimensional (2D), with drift 
vortex motion taking place both in the directions 
parallel and perpendicular to FL{t). The ratchet 
effect is one-dimensional (ID) only if ^ = 0^, 90^ 
, that is, when ^^{t) oscillates along directions 
of high symmetry of the rectangular dipole lat- 
tice {x and y directions, respectively). In these 
cases drift vortex motion takes place only in the 
direction of the driving force. 

Next the results of the numerical solution of 
Eq. are reported for the following parameter 
values: d = zq = 2(^, X = 10^, ax = 20^, ay = 40^, 
m = 0.25^0^0- 

3.1. One-dimensional ratchet 

Here the ID ratchet effect is considered for = 
0, that is Fi,(t) oscillating in the x— direction. In 
this case the vortex drifts only in the x— direction 
{Vy = 0), and Vx depends on Jt, and P. Typ- 
ical Vx vs. Jt curves for < a < 90^ are shown 
in Fig. IHl In this case the vortex drifts in the 
negative x-direction. For 90^ < a < 180^, the 
vortex drifts in the positive x-direction, and the 
Vx vs. Jt curves are identical to those for 180^ — a 
with Vx replaced by —Vx- For a = 90^ there is 
no drift vortex motion, and Vx = 0. The most 
important characteristics of the Vx vs. Jt curves 
for < a < 90^ are the following. Drift vortex 
motion only occurs if Jt is larger than a mini- 
mum value, denoted by Jo(q^). For Jt > Jo{(^) 
the Vx ys. Jt curves have two distinct regions: 
one for Jo{a) < Jt < Jm{<^)^ and another for 



Jt > Jm{oi)^ where Jm(<^) denotes the value of 
Jt for which Vx is minimum (| | is maximum). 
It is found that in the region Jo(q^) < Jt < Jm{(^) 
the vortex displacement during each period is an 
integer multiple of the lattice parameter ax. The 
Vx vs. Jt curve consists of a series of plateaux 
where Vx is independent of Jt. In each plateau 
Vx = iax/P,i = 1,2,3,.... As shown in Fig. 
the plateaus are clearly visible in the curve for 
P = 1200r. For P = 24000r the plateaux are 
too narrow to be distinguished in the scale of the 
figure. For Jt > Jm{(^)^ the Vx vs. Jt curve is 
strongly affected by the time dependence of the 
driving force. For the square wave, Vx drops to 
zero very quickly with Jt, while for the sinus wave 
Vx drops to zero slowly and is a complicated func- 
tion of Jt. The dependence of Vx on a for fixed 
Jt is shown in the inset of Fig. [Hlb, and will be 
discussed in more detail later. It is found that for 
both time dependencies Jo{a) :^ Jc( 180^ — a) and 
Jm{<^) — Jc{o^)- More details of the dependence 
of 14; on Jt, a and P is given next. 

a = 0: In the steady state regime there is no 
motion in the y— direction. The vortex moves 
only in negative x-direction, along a line connect- 
ing the minima of Uy-dai'^)^ for instance y = ay/ 2 
(see Fig. inia). Thus the vortex dynamics is one- 
dimensional in the potential Uy-da{x,y = CLy/2) 
for a = 0. As shown in Fig. this ID po- 
tential has characteristics similar to ratchet po- 
tentials used in other contexts Q]. For Jo(0) < 
Jt < Jm(0) the vortex moves only during the 
half periods of oscillation when FL{t) points in 
the negative x-direction. When FL{t) points in 
the positive x-direction the vortex is pinned be- 
cause Jm(0) ^ Jc(0) (see Fig.m. For Jt > Jm(0) 
the vortex moves both in the negative and posi- 
tive X— directions, but the displacement is larger 
in the negative x— direction due to the lack of in- 
version symmetry ofUy-da{x^y = ay/2). 
< a < 90^: Vortex motion is found to be 
two-dimensional, but drift takes place only along 
the negative x— direction. In the ^-direction the 
vortex oscillates with zero average. Drift vor- 
tex motion is similar to that for a = 0. For 
Jo (a) < Jt < Jm(<^) the vortex is pinned when 
Fl(^) points in the positive x-direction because 
Jm{o^) — Jc{o^). However, since Jc(q^) decreases 
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Figure 5. Vortex drift velocity (in units of eo/r?) 
vs. Jt for driving force parallel to the x-axis {0 = 
0), and magnetic moments at an angle a with the 
X-axis. Labels: SIN and SW = sinus-wave and 
square- wave driving forces, respectively. Inset in 
b) 14 vs. a for Jt = 0.5 and P = 1200. P in 
units of r = r]^/eo. 



as a increases, the value of Jt for which Vx is min- 
imum ( Jt = Jm{(^) ) decreases with a, and the 
value of Vx at the minimum increases (| 1^ | de- 
creases) with a. As shown in the inset of Fig.IHlb, 
for Jt fixed, the value of Vx is weakly dependent 
on a as long as Jt < Jm{(^)' For Jt > Jm{(^)^ 
Vx drops to zero rapidly with a. This result 
can be understood by noting first that the av- 
erage slopes of the Vx vs. Jt curves in region 
Jo (a) < Jt < Jm{o^) are essentially independent 
of P, as shown in Fig. Second, the average 
slope of the Vx vs. Jt curve is linear, except for 
Jt close to Jo(q^), and weakly dependent on a. 
This occur because the driving force is large com- 
pared with the pinning force in the direction of 




-25 X 25 



Figure 6. One-dimensional vortex pinning po- 
tential Uv-da{x^y = %/2) (in units of eo<i), for 
a = 0. X in units of ^ 



drift. An estimate for Vx in this case is obtained 
by neglecting the pinning force during the half pe- 
riod of oscillation where vortex drift motion takes 
place. The result is Vx ~ —0.5JT/Jd{^o/v) 
the sinus wave, and Vx ~ —0.8JT/Jd{^o/v) 
the square wave. These estimates are in reason- 
able agreement with the values of Vx shown in the 
inset of Fig. [Sib for Jt < Jm{o^) 

Tuning the ID ratchet effect is illustrated in 
the inset of Fig. ^h. Basically, by rotating the 
magnetic moments the vortex drift velocity can 
be switched off, on, and reversed. 

3.2. Two-dimensional ratchet 

Now the 2D ratchet effect for ^ 7^ is consid- 
ered. In this case V depends both on and a, 
as well as on Jt, in a complicated way. Typical 
results for the dependence of V on Jt and 6 is 
shown in Fig. EJ and Fig. [HI for the square- wave 
driving force. The curves for Vx and Vy vs. Jt 
are similar to those obtained for ^ = 0. There 
is a value Jt = Jo below which Vx = Vy = 0^ 
and a value Jt = Jm for which both Vx and 
Vy are minimum. It is found that Jm — Jc(| 
9 - a I), and Jo ^ Jc(180^- \ 9 - a \). For 
Jq < Jt < Jmi the Vx vs. Jt and Vy vs. Jt 
curves show plateaux where both Vx and Vy are 
independent Jt- Each plateaux correspond to a 
vortex displacement by a lattice translation vec- 
tor, R = Uxaxyi + riyayf, Ux, riy = 0, ±1, ±2, .... 
The direction of V (angle 7 with the x— axis) is 
not, in general, parallel to the direction of driv- 
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Figure 7. Vortex drift velocity components (in 
units of eo/r]) vs. Jt for magnetic moments par- 
allel to the X-axis (a = 0), and square- wave driv- 
ing force at angle with the x-axis and period 
P = 1200r. Insets: angle between the vortex 
drift velocity and x-axis, 7, vs. Jt- Horizontal 
lines: 7 = 6> + 180^. 



Figure 8. Vortex drift velocity components (in 
units of eo/f]) vs. Jt for magnetic moments at 
a = 30^ with the x-axis, and square-wave driv- 
ing force at angle 6 with the x-axis and period 
P = 1200r. Insets: angle between the vortex 
drift velocity and x-axis, 7, vs. Jt- Horizontal 
lines: 7 = 6> + 180^. 



ing force oscillation, as shown in the insets in 
Fig. and Fig. [HI However, as shown in Fig. 
□a ( a = 0, 6> = 11.25^) and Fig. [Hla. {a = 

30^, = 15^), there are regions where the two 
directions nearly coincide (7 = 180^ + 0). This 
occurs because for Jt in these regions the driv- 
ing force is large compared to the pinning force 
during the drift vortex motion. Similarly to the 
case ^ = discussed above, the average slopes of 
the Vx and Vy vs. Jt curves are linear in these 
regions. For a = 0, = 33.75^, 7 differs consid- 
erably from that of driving force oscillation. For 
a = 30^, 6 = 30^, there is also a region where V 
is parallel to the direction of drive, but the aver- 
age slopes of the Vx vs. Jt and Vy vs. Jt curves 
are not linear. In this case the driving force os- 



cillates in a direction parallel to that of the mag- 
netic moments, but the dependence of the drift 
velocity on Jt differs considerably from that for 
a = 6 = 0. This shows that V depends both on 
6 and a. Tuning of the vortex drift velocity by 
rotating the direction of the magnetic moments 
is also possible for ^ 7^ 0. The dependence of Vx 
on a for fixed Jt is, essentially, similar to that 
shown in the inset of Fig.[5Jb, since the Vx vs. Jt 
curves for ^ 7^ have the same basic structure as 
that for 6 = 0. 

4. Discussion 

First, an order of magnitude estimate of the 
vortex drift velocities obtained above is given. 
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Consider the case ^ = and assume that ^ = 
20nm, and r] = r]od, with rjo = 7 x 10~^ Nsm~^ 
12]. In this case d = zq = 2^ = 40nm, A = 
10^ 200nm and m = O.250o^o ^ 10 
(/iB =Bohr magneton). According to Sec. Ol 
the maximum values of | 14 | are | Vx \max^ 
0.bJc{a)/ Jd{eo/r]) ~ 10Jc{a)/ Jdm/s for the si- 
nus wave and | \max^ 0.8Jc{a)/ Jd{eo/r]) 
16Jc{a)/ JdTn/s for the square wave. This gives 
for the sinus wave | Vx \max^ 10m/ s for a = 0, 
and I Vx \max^ "Irn/s for a = 45'', and | Vx \max^ 
16m/ s for a = 0, and | Vx \max^ 6m/ s for 
a = 45^, for the square wave. These velocity val- 
ues are of the same order of magnitude as those 
reported in Ref. • The value of the unit time is 
r = r](,/eo ~ 10~^5. Thus, the transport cur- 
rent oscillation frequencies, u = 1/P, used in 
the calculations are u = 800 KHz, 40 KHz for 
P = 1200r, 24000r, respectively. 

The results obtained in this paper are be- 
lieved to be representative of low- Tc supercon- 
ducting films with magnetic dipole arrays placed 
on top. First, the particular set of parameters 
used, (i ~ ~ ^, are typical ones. For instance, 
in the experiments with arrays of magnetic dots 
with permanent magnetization placed on top of 
superconducting Nb films, reported in Ref.|15j. 
d = 20nm ~ ^. The magnetic dots are separated 
from the film by a thin protective layer of thick- 
ness ~ 20nm, so that the distance from the mag- 
netic dipole to the film is zq ~ ^. Second, since 
the vortex drift velocity depends on the scaled 
parameters m/^o^o, d/zo, J/Jd^ many supercon- 
ducting film-dipole array systems are equivalent. 
Third, the results obtained in this paper do not 
depend on the particular dipole array used. What 
is essential for the ratchet effect is the periodic 
vortex pinning potential with broken spatial in- 
version symmetry created by it. Most dipole ar- 
rays will do, since the interaction between the vor- 
tex and a single dipole lacks inversion symmetry. 

The London limit is valid for vortices in low- Tc 
films. However, when a magnetic dipole is placed 
close to the film, it certainly breaks down if the 
dipole field destroys superconductivity locally in 
the film. Roughly speaking, London theory is 
valid as long as the maximum dipole field at the 
film is less than the upper critical field, that is. 



m/z^ < ^o/(27r^^), or m/(^o^o) < (^^o/OV^tt. 
For the parameters used in the above calculations 
{zq = 2^) this gives m/(^o^o) < 0.64, which is 
larger than the value of m used in this paper. 
The London limit would be a better approxima- 
tion if the present calculations were carried out 
for larger values of zo/(_. However, the results 
would be identical to those described above if m 
and d were scaled by the same factor as ^o/^- For 
instance, if zq 2zo, Jc/Jd would remain the 
same \i d ^ 2d and m 2m, but the upper 
limit of m/(j)oZo for the validity of the London 
approximation would increase by a factor of 4. 
The present model also breaks down if m is suf- 
ficiently large to create vortices in the film. The 
threshold value of m for spontaneous vortex cre- 
ation, estimated as m ~ O.T^o^o using the results 
of Ref. 221: is larger than m used here. 

The results obtained in this paper for a sin- 
gle vortex also apply to vortex lattices pinned by 
dipole arrays if the vortex lattice is commensu- 
rate with the dipole array and if the vortex den- 
sity is at most one vortex per dipole. In this case 
the driving force is expected to move the vortex 
lattice as a whole, that is preserving the spatial 
order, so that the effects of vortex- vortex interac- 
tions are negligible. 

In conclusion then, this paper demonstrates the 
ratchet effect for vortices in a thin superconduct- 
ing film pinned by periodic arrays of magnetic 
dipoles placed on top of the film, and shows that 
the ratchet effect can be tuned by rotating the 
magnetic moments. In the case of magnetic mo- 
ments free to rotate, the ratchet is created and 
tuned by a magnetic film applied parallel to the 
film surfaces. 
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